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SHOCK  WAVES  IN  SOLIDS 
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(Moscow  Higher  Technical  School  imenx 
K.Ya.  Bauman }  submitted  to  editors 
16  April  1958) 


.  •  ^  f  <r'Hr*iric  waves  in  solids  under 

The  propagation  of  shock  **>-v  e-  ..  .H+v  _f  the 

high  pressures  co»siderinS  the  co«pre..i  ^  .^on  of 

is  exasBxned  xb.  tJriXe»  articxe- 

“e  miteri.l  is  studied  in  the  so»«  of  “tnordent shock 

wave  8 rsd  the  reflected  wave  zone,  wnei  e  the  body 

finite  dimensions .  The  solution  is  obtained  x«  th« 

f0r*  iotio“rof  t?.;  incident  .hoe*  •»  »  investi¬ 

gated  by  means  of  «  special  soiution  of  differential 
equations  of  motion,  turn  wave  ‘  f  dif„ 

surface  is  investigated  by  a  general  soln.ja  of  dx 
f erential  equations  of  -tion.  The  results  obtained 
make  it  possible  to  determine  all  of  ^e  pax ameter 

b0tiliS^^t^i«"fthfthJrIS«l  results  obtained 
.  dynamic  W  deformation  stress  can  be  derived 

with  the  appropriate  experiments. 


To  solve  the  dynamic  problems  arising  when 
detonation  or  shock  charges  are  applied  to  a  solid 
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it  Is  necessary  to  know  the  dynamic  deformation— stress 
curve*  An  experimental  determination  of ^ such  diag- 
t a s  ' pr 0 b o ti is  c o n $ i d. ® r a b JL 0  d  1. f f  1  c u .1 1 ie$ #  X  t  is  kilo wb 
t'h^t  when  a.u  instantaneous  pressure  kis  applied  to  a 
solid, -which  then  falls  in  the  course  of  time t  elastic , 
plastic  or  shock  wave*  can  be  propagated  in  the  solid* 
depending  on  the  Esagn.it vide  of  the  applied  pressure  ami 
t he • properties  © £  t h e  r © 1 i d  i & e e  wo rk  ( 1 ))# 

In  the  present  work  the  solution  of  the  dynamic 
problem  at  high  detonation  charges  is  analyzed*  where 
shock  waves  are-  propagated  in  a  solid  medium*  In  this 
a  r  e  a  of  gr  eat  dyn  am  i  c  p  r  e  a  s  u.r  ©  s  «  wh  e  r  e  t  h  e  a  olid  modi  urn 
behaves  like  a  compressible  fluid  (1)  *  the  dynamic 
deformation-stress  volume  diagram  is  only  . qualitatively 


known 


it  can  be  given,  for  example*  by  the  egua 


tion  having  the  .foil. owing  form*  6‘  ~  ffa  «*  —  Sc 

where  A. *  6*  *  are  unknown  constants *  •  In  such  an 
approximation  of  the  S'  ~  g  curve*  one  can  obtain  in 
the  form  of  finite  formulae  the  laws  of  velocity  and 
pressure  change  in  hath  the  ir.etckmi  shock  wave  and  the 
r e f 1 ecte d  wa ?e  { we r t h 0  b o dy *  t o r  i msi an ce,  to  be  a 
plate  of  finite  thickness) *  St  is  known  that  if  impact 
loading  is  applied  at  one  side  of  a  finite  body  (a 
plate)?  at  the  rear  a  part  of  the  material  will  be  split 
off  with  great  velocity*  Using  this  phenomenon*  it  is 
possible  by  means  of  the ' solution  obtained  in  this  work 
to  find  the  unknown  constants  A  *  %  and  *  the  dynamic 
da  forma t i on- a tr a a 3  volume  curve  with  compression  will 
be  determined  by  this* 

We  will  examine  the  propagation  of  a  shock  wav# 


in  a  solid  *  w 

here  a 

great  i: 

n  s  t  a  n  t  a  n  e  cm  s  p  r  e  s  0  nr  e 

which 

$  ufe  s  e  quen.t  X  y 

do ore as 

as  with. 

tim®  Ik  applied  to  the  let- 

ter  * 

is 

knO'rOi.  t 

hat  ths 

n  do  forma  t  i  on-  sires® tf 

com- 

pr  e 3 s i on  0 urv 

e  for  s 

oli&B  h 

the  form  shown  1b. 

Fi&»  1 

( iii©  work  Cl) 

> 

/ * 

Pig*  1 


If  the  pressure  in  the  wave  exceeds  the  pressure 
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corresponding  to  the  discontinuity  b  of  the  6  -  6  curve, 
shock  waves  will  be  propagated  in  the  solid.  It  is 
apparent  that  under  these  pressures  the  solid  may  be 
regarded  as  a  quasifluid  (1), 

It  is  known  that  the  basic  equations  of  one ^dimen¬ 
sional  motion  have  the  form 

da  ^  dp  _  ^  da  _  dv  ^  ^  ^ 

at  1  dh  '  dh  di 

where  u  -  is  the  particle  velocity;  v 
volume 5  t  —  is  time;  p  -  pressure;  h 
coordinate  * 

If  one  introduces  the  terms  tS 

-  1, 


we  obtain 


du  _  d<s  ^  da  ch  <iu  fh  '  f 

dt  dh  dz  dh  dh  °  di 

The  equation  of  state  (curve  fe  ~  e  1b.  Fig*  1)  is 

approximated  by  the  relation  having  the  following  forint 

ci  ~  0<D  rr.  A  (e  —  ea)  ~3t  (  4  ) 


*  tl 

9 

then 


is  specific 
is  the  Lagrangian 


P  t 


6 


(2) 


where  ,  A  and  c&  are  constants* 
¥  e  then  o h  t a i n 

~de  (*  — «o)+4  ’ 


(5) 


The  velocity  of  sound  in  this  case  is  determined 
by  the  expression 


(«+ 1)  i/.i 


£  4- 


Po  d  £  (€  —  £$)* 

V  3 A  po _ 

~  (l  +  «o) 


■i/^C  - 

*  Po 


(6) 
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The  velocity  of  propagation  of  disturbances  from 
particle  to  particle  will  fee  equal  to 


We  write  equation  (3)  1b  the  fora 

du  _  d  0  d  t  3 A  0{s— &$) 

dt  d  $  Bh  (t  — e&  !4  dh ' 

_  .. 


(8) 


We  will  find  particular  solutions  for  this  system 
of  equations*  We  will  assume  for  this  purpose,  as  usual , 
that  g  =  6  (u)  9  then  we  will  arrive  at  the  aquations 


d u  __  d s  ^ if  ^ 

w~~~inr  ~Tk  '  ■  cn  du  St  *  ’ 


whence  we  obtain 


(9) 

(10) 


The  solution  of  equation  (10)  has  the  form 
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may  fee 


tera  of 
lad&p&S 


th&n  6 
form 


h  =  \f"%  ------  f  +  F  M  =~  P*  **  +  W- 


(M) 


ff  * 


In  the  equation  of  state  (4)  the  specific  solution 
written  in  this  form 


A «  4- 1 


i..  i  /  34  . — 1 - 1-  const 

-  y  -j--  ■»-««,■ 

/  3,4  _  .4_.  F(ft) 


4-  i/  ..if ’JL_  -f-  ; 

'W  |r  (*  —  £*)* 


(12) 


We  will  now  seek  the  g^a^ral  sv>£u  lions  of  the  sys* 
equations  ( 8)  *  Having  ;La?e rted  the  dependent  and 
3  dent  variables  *  we  arrive  at  the  equations  j 
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#  A 


6?  e 


2r  / 

4  J*L 
d«r 


$  t 


e^,w  *  'rJ  '  TfHKnnftMMPl^ 

*  *  ib 

a  a 


(13) 


We  find  from  this  that 


(P  ti  ____  #*/  ,  d*/ 

“  dsda  d>*1  d# 


In  the  equation  of  state  (4)  w@  will  have 


#/_  ^  _ 3.4 

d(st  —  sa)s  pefa  —  ecy  &M* 


(!5> 


We  will  introduce  the  new  variable  e  «  -g%~ 
-  and  equation  (15)  assume*  the 


f L + JL  d?  __  34_  a8 * 

dz*  2-  d  z  fca7  da* 


{ 16) 


Assuming  that  a 


2  ^  3A 

p  n  ,  then 


5 


rjr  _ 

Po  E 


e  - —  $ 


:  '■ —  ®  (e  —  eo)  ~ 

e  -{*  | 


a  or 


w 


2 


a 


equation  (16)  assumes  the  form 


d2(zf) 

dz2 


9Hz  t) 

d  az 


The  solution  of  this  equation: 

i  ~  'Ml 


",d  P2  arbitrary  function, 

according  to  ti 


IV e  land  tram  the  2nd  equation  (13}  that 


f  r 
~dZ 


z7  dt 
a  &  z 


-j-  f  £i  ,F 3  ..1  4  ^2  ]  J'Z  __  _| 

z  ^  ~2  J  <?  ft 


whence  it  follows  that 


A 


fjy 

a 


{ 2 


?.)  -f  0~\  +• t\)  ■ 


If  the  function  d>  -  m  ( c  „\  ■ 
satl5fi„  the  second  eolatiS.  (13)  |  then""* 


d<4' 


dy 


* ■  ^4* 

n  p0  ~r~' 

d  t 


Pc* 

a 


djP 

dz 


6 


(17) 

m ) 

(tha 
*  argu- 

d  h 

~dV  ' 

09) 

which 

OT) 


J 


rhere 


Z 


rhere 


Whereupon  the  solution  of  <l8)  assumes  the  form 


fl  and  f2  are'  nevr  a^itrary  functions, 


r  + 

fi  +•  h 

- 

r 


p» 

a 


.  <22) 

. . ..  ..  J 


ihe  special  solution  of  (11)  assumes  the  form 


&*<* 

a  d  z 


fW  4® 

W^S>/ 

du 


~  T  ^(«)  . 


(-3) 


Whereupon 

d  u  ■ 


r- 


J?  c 

Pc  dz 


dz„ 


(24) 


We  multiply  (2p)  by  du,  then 


f fj/"- -*-  —  +  *1/ 


jvt_  Jda  _____ 

j»  8  z  3  w 


fife  —  F/e/)  du, 


Whence 


?/ 

r 


it 

fl 


2 L</* 

dr 


^  -  ii )  d  u  or 


~  i 


*  §  (a)  du,  (25) 


7 


rtzn 


Hsncfet  along-  the 
which  are  determined  by 


efa&ra.ctiyrlstz cs  in  plane  (u$ 
the  b p  Si  c i  a  1  s o-l  u t  i on 


«) 


(26). 


If  W(n)  3«  0 v  then'  4*  *  const,  >■  since  without  limit 
ing  the  general  representation  it  may  always  he  assumed 
that  ecm&t  ♦  a:0 *  hence  if  P(xx)ar O,  along  the  characteristic 
di  «s  0  (2-7)*  In  the  equation,  of  state  (4)  *  we  have 


tt  =s  *ir  x  4“  const* 

il:>  £  4  -  h  f  f } .  ( 2  8 ) 


ft 


,x..  CJL*. 


■L  f 


Fft)  z 


If f  for  example*  we  deter Mine  ib  along  the  charac¬ 
teristic  "a  *  z  s«  const  *t  then  in  this  individual  case, 
where  F(sk)  a?  0*  £  a  0* 

If  *|J  is  determined  along  the  characteristic 

v<  ~  a  as  const.,  then  when.  F{»)  »  0,  t"  «  0* 

l 

We  will  look  at  the  solution  of  a  concrete  example 
Let  us  assume  at  the  moment  of  time  t  ~  0  in  cross 
section  h  s  0  there  is  a  pressure  p  -  pH  applied*  where 
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pKis  greater  than  6%  ,  where  6&  is  the  pressure  corres¬ 
ponding  to  the  discontinuity  point  in  curve  6  -  6  (Fig.l) 
which  then  is  found  according  to  the  law 


p-p«^{  i _ Yn 

Ph  -  Pa  l  ) 


(29) 


where  p# , x  and  n  are  constants. 

It  is  known  from  detonation  theory. that  it  is 
possible  to  assume  p0  a  0,  n  »  I  and  X  =  s  where  1 

is  the  height  of  the  loading  and  D  is  the  velocity  of 
the  detonation. 

We  find  the  law  of  Motion  of  the  medium  in  the 
passing  wave.  The  motion  of  the  medium  in  the  wave  will 
be  subject  to  the  specific  solution  (28)  where  it  is 
'necessary  to  determine  F(a): 

u  —  z  .-const 

h~  F(z) 

a 


We  v/ill  determine  the  constant  in  the  first  equa¬ 
tion  of  the  system  (A) * 

When  n  »  \in  ,  z'  = s  zH  ¥  inasmuch  as  at  the  wave 

front 


form 


Using  equation  (4)  , 


{29)  can  he  written  in  the 


(30) 


The  first  equation  (A)  now  assumes  the  form:  u  ~  z  -  zQ  • 

(31) 


For  the  initial  moment 


9 


<32 


« 


«  - *-«  *0- 


If  it  is  assumed  in  equation  (A)  that  h  «  0,  and 
ie  time  t  is  expressed  from  (30) %  then  we  determine 


(*)s 


m 


H  **  l 
a 


l 


A  z 


f- 


(33) 


System  (A)  is  finally  written  in  the  form 

u  =  a  *  s«  „  , 


hi  \  *  /  J; 


(34) 


If  n 

Lve  to  55  5 


I ,  we  derive  the  quadrat it  equation  rela 


T 


« -r 


ah 

fi0(f  ~H) 


0. 


The  solution  of  this  equation  will  be 


Z  s=  tf-HV 


z„t  +  |/'J^2t*4-4%h  A((  “FA 


(35) 


We  find  the  law  of  motion  of  the  shock  front 
h(t)  : 


IX 


dh  u  ~jr  c  ~f_fo  z  —  z0  ~4~ 

—  ~~  -  -  2 


(36) 


sere 


D  is  the  velocity  of  the  shock  wave , 


c  is  the  velocity  of  sound  in  metal. 


10 


We  express  the  local  velocity  of  sound  c  by  at: 


(** 


±J2 


:'i  -j~  a 


a 


■  u.  ■» 


Thus »  the  velocity  of  wave  B  will  be  equal  to 


/)v  =- 1 4 


C0  “•  .  %  4 1  Tg 

2  •  ‘  2a 


(38) 


By 

means  of 

equation  (35)  ««  eliminate 

z 

...  c0  —  r0 .  |  « 

Pc  ^ 

2  4  t) 

}-[i  4 

gp-f-'i 

2h  2  1 4 rw*  t*  4-  4  v0ah  (t  -I- 

ga 

t  X  j  2{i  4pt} 

(39) 


In  t  e  gr  «  i  ±el  g 
we  find  1;ia«2  lav  of 


this  equation  which  when  h  «  0, 
propagation  of  the  shock  wave  s 


f 


0, 


C40)- 


If  the  body  has  finite  thickness  d ,  then  when,  the 
wave  reaches  the  t&ck  side  of  the  solid,  ix  will  be  ref  lec 
ted  from  the  free  surface  as  a  rarefaction  (expansion 
/"“Russian:  discharge^)  wavs,  moving  toward  the  shock 
wave.  The  motion  of  ""the  medium  in  the  rarefaction  wave 
is  determined  fey  the  general  solution  (<Js)  •  The  arbit¬ 
rary  functions  in  the  general  solution,  are  determined 
primarily  from  the  condition  cf  the  connection  between, 
the  reflected  and  incident  wave® 5  secondarily,  by  the 
fact  that  When  t £T;  h  -  p6d,  •T  «  0,  where  d  is  the 
thickness  of-  the  plate.  The  first  condition  implies 
that  on  the  lino  «.  ss  *  a*  takes  place  according  to  (25) 
and  ( 33) - 


|  4>{-;>)du 


4  frY4 

I 


H 


2-t 


F 


i  I  £& 


!  -f- 


« — 1 


/  y,lh" 

/  JL« 


11 


r 

If  n  ~  1 ,  the»  — 1  -f* 


!n 


or 

—  t  2'2  -f-  xzHZ  Jn  ----r  -- /,  -f  fi , 

*r 


f41) 


where 


f i  -  (2sb  -  35*)  and  f 2  ~  f?(-z0). 


x\  Is  the  integration  constant* 

The  second  condition  provides  that  when  h  s  pQd 
the  condition 


ad  -~ji  ■{-  ft  Jr  i  fz  —fj 


is  fulfilled «  where 


Thus  f  by  Keans  of  {4.1}  and  (42)  the  parameters  of 
the  reflected  wave  can.  be  determined*  The  found  solution 
is  complicated,  for  investigation ,  hence  we  will  find  a 
less  precise,  although  considerably  simpler  solution  of 
the  problem*  It  is  known  from  detonation  theory  that 
the  pressure  on  the  surface  of  the  plate  falls  according 
to  the  law 


where 


(43) 

(44) 


The  velocity  of  the  shock  wave  D  is  determined 

y 

by  the  formula 


12 


V  * 


u  “4*  £  *4~  £q _ _ 


We  take  the  equation  of  state  in  the  form 


p  ~~  A-t/  —  8< 


where  and  k  are  constants, 

The  motion  of  the  shock  wave  will  be  determined  fey  special 
solutions  of  the  system  of  equations  (1)  written  in 
Eulers  form: 

„  „  ,  h — I 

r— —  a.  (47) 


Jfcssfii-f  Ff« j. 


The  velocity  of  the  shock  wave  By  is  cletermixied 
by  the  expression 


D¥  3=: 


Function  F(u)  is  determined  by  the  condition  that 
when  x  s=  0  the  pressure  varies  according  to  the  law  (43) 


P  ~~ 


,i_\* 
ot  j 


whence  follows  that 


•£L  =  /  J&lVa  ss  Ja 

/  U  / 


where  is  the  initial  speed  of  sound  in  the  plate 

where  p  »  p  s 

f« = /If  =ir**r(&±*p  .  (5 

By  means  of  (47)  equation  (50)  assumes  the  form 


15 


I  l  «  \ 


2 


D 


'52) 


When  x  s  0  tre  determine  Hu)  according  to  (48) 
and  (52); 

v  /  <  *“H  .  \f  ,  *~1  \~‘  / 

I  2  A  2  /  Z>c* 


Equation  (48)  now  .has  the  form. 


X' 


Jfe  + 1 


f  — L 

*■  d  i 


k — i  \ 


Q?  i — ~~~  w 

*  / 


/ 

j 


(53) 


When  k  -•=  3  this  expression  is  written  in  Ihe  form. 


/ 

x  4-  2 a )  (  f 


<4  „Lj 

/i  ~  &  / 


(54) 


We  convert  this  expression  to  the  form 


X  (u  4-  c„)  -f-  2i/) 


Dr 

i 


•  («  -f  4)  —  <4  ]  ,  (55) 


or 


0  /j*  ,  .  r  w  5  d 

2  - <r  4-  //  - 3c,,  ~  x  — 

/  [  /  "  l 


2c,  |  4 


1  .  Dt 


cl 


C H  C &  **'  X<  , 


D 


0. 


2 

Inasmuch  as  u <£c  ,  neglecting  the  term  witn  u  , 

o 


we  find  that 
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p 

f 


H  j 


X 


Dt 


u  = 


Dt 


v>  “ 


x  D  n  cH 


*  a 


When  k  «  3  equation  (4*9)  assumes  tho  farm 


r,  dx  .  ' 

Lh>  ~ - U  ~i‘  ca 


ir 

L'fi  ^ 


Dt 


cK  1 


{ . . £aj 

X  'll.,}  V*v  iT  ^ 


(S7) 


We  assume 


n# 


“  1 4-. 

2c~ 


£  _ 

v  - — 


2D  l 


we  then  get 


£1 

dx 


ik. 

6 


2": 


D 


(58) 


We  will  designate  the  reJaUon 


2: ,  then 


or 


dz 

d  x 


■.  r. 


it 

15 


d  m  x~ 


J> 

o  — 

O 

-*■’*■*  / 

c* 

Dz 

(  3 

\ 

D 

'C9 

j  dz 

ZZL 

_ _ _ _ . _ ^  ■*r 

£ZZ 

F>y  l 

~fj~~ 

... 

5  ’  l 

~  i>" 

dZ  -4- 

1  a 

z* 


At  i  e  r  i  n  t  e  e  r  a  t  i  o  n  we  o  b  t  a  I  n 


the  law  of 


Front  this  we  determine  x  &  x(t)  - 
motion  of  the  shook  wave- s 


where 


X  XSA 


±  %/  il,  Vttn 

r  4  \7) 


B 


r?- 

V/j 


if!"'  \  $ 

*  *  a  \  *  /?  ? 


£«_ 

Li 


1ca  \ 


\  /* 


The  obtained  relations  cast  be  need  to  determine 

pressure  in  the  parsing  waves* 

To  determine  the  pressure  in  the  passing  wave  at 

a  given  distance  x  «  h  one  mus fc  determine  according 

to  formula.  (6-3)  when  x  a  h9  then  according  to  (56)  deter- 

mine  tr  4  in.  accordance  with  which  o,  is  determined  with 
n  h 

(V?)*  The  pressure  is  determined  by  iweans  of  the  e  qua  - 

t i on  0 £  $ t  a  t  e  { 4  6  !  * 
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